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Abstract. In this paper we have constructed an approximation for the 
Harris flow and the Arratia flow using a sequence of independent stationary 
Gaussian processes as a perturbation. We have established what should be 
the relationship between the step of approximation and smoothness of the 
covariance of the perturbing processes in order to have convergence of the 
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approximating functions to the Arratia flow. 
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^ ; Introduction. It is well-known |TJ, that the solution to the Cauchy 

problem for SDE 



| dx{t) = a(x(t))dt + b(x(t))dw(t) 

\x(0) = Uq 



> 

O 

with continuously differentiable functions a, b having bounded derivatives, 
^ ■ can be obtained via discrete time approximation. Namely, if we define the 

sequence {x™} by the rule: 



X™ = X £l, C+l = C + -«) + A=K*n)Zn (2) 

m \/m 



where n > 1} is a sequence of independent standard Gaussian random 
variables, then the random functions 



x m (t) = ml — --t)xT+m( t )xT+i, t G 

V m I \ ml 



k k + 1 



m m 



, k = 0, . . . , m—1 



weakly converge in C([0, 1]) to the solution of (GQ). 

In this paper we study similar to (ED difference approximation for coa- 
lescing stochastic flows. As is known [2], such flows are not generated by a 



Gaussian white noise in the space of vector fields. In order to understand 
how the flow with coalescence is arranged we can consider its difference 
approximation. As a perturbation we select a sequence of Gaussian sta- 
tionary processes. In order to allow the coalescence of the trajectories of 
individual particles in the limit, the covariance functions of these processes 
are chosen to be less and less smooth at the origin. On the other hand, in 
order the limit flow to preserve the order, the step of approximation must 
be sufficiently small. The relationship between the step of approximation 
and smoothness of the covariance of the perturbing processes explains to 
some extent the structure of singular stochastic flows. 

1. SDE and stochastic flows on the real line. The main object 
of the article is the Harris flow of Brownian motions on R. Let <p be a 
continuous real positive definite function on R such that (p(0) = 1 and ip 
is Lipschitz outside any neighborhood of zero. 

Definition 1. The Harris flow with cp being its local characteristic is a 
family {x(u, •); u G R} of Brownian martingales with respect to the joint 
filtration such, that 

1) for every u\ < and t > 



It is known that the Harris flow exists [3]. If the function ip is smooth 
enough, the Harris flow can be obtained as a flow of solutions to SDE. 
Namely, for a sequence of standard Wiener processes {wt] k > 1} consider 
the following SDE 



x(ui,t) < x{u2,t) 



2) the joint characteristics are: 



d(x(ui, -),x(u2, -)){t) = <p(x(ui, t) — x(u2, t))dt. 



00 




(3) 



k=l 



where a = (a,k)k>i is a Lipschitz mapping from R to I2 such that 



00 




2 



and 

oo 

^2a k (u)a k (v) = ip(u - v). 

k=l 

Then the flow corresponding to (ED is the Harris flow with the local 
characteristic </?, and furthermore it is a flow of homeomorphisms. Note, 
that the Harris flow could be coalescent [3j and, in this case may not be 
generated by SDE. By this reason it is interesting to consider discrete 
approximations for the flow built in a similar way as approximations to 
SDE. Consider a sequence of independent stationary Gaussian processes 
{£,n{u)' } iiGl,n> 1} with zero mean and covariation function T. Suppose, 
that T is continuous. Define a sequence of random mappings {x n ; n > 0} 
by the rule 

x (u) = u, x n+ i(u) = x n (u) + £ n+1 (x n (u)), u e E. (4) 

Note, that the continuity of T implies that the processes n > 1} have 
measurable modifications. This allows to substitute x n into £ n +i- The 
independence of > 1} guarantees that £ n+ i(x n (u)) does not depend 

on the choice of these modifications. We will need the following description 
of one and two-point motions of {x n ; n > 0}. 

Lemma 1. The sequences {x n (u);n > 0} and {x n (u2) — x n (ui);n > 0} 
have the same distributions as the sequences {y n (u); n > 0}, {z n (u); n > 0} ; 
which are defined by the following rules: 

Vo = u, Vn+l =Vn + rjm 

z = u 2 - iti, z n+ i = z n + y / 2r(0) - 2F(z n )r] nj 
where {r] n ; n> 1} is a sequence of independent standard normal variables. 

The proof of the lemma can be obtained easily by calculating conditional 
distributions of x n+ \ under given xo, . . . , x n , and is omitted. 

It follows from Lemmadl that the sequence of random mappings {x n ; n > 
0} is similar to the Harris flow. All its one-point motions are Gaussian sym- 
metric random walks. But the mappings x n for n > 1 are not monotone. 
In the next section we will prove that any m-point motion of {x n ; n > 0} 
approximates the m-point motion of the Harris flow. 
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2. m-point motions. In this section we will consider the limit behavior 
of x n under a suitable normalization. Let us define the random functions 

x n (u, t) = n (—— ~ tj Xk(u) ^ n [ t ~~) x k+i(u), 



«e1, te 



k + 



, k = 0, . . . , n — 1. 



n n 

Our first result is related to the n-point motions of x n . 

Theorem 1. Let T be continuous positive definite function on R such that 
T(O) = 1 and T has two continuous bounded derivatives. Suppose that x n 
is built upon a sequence k > 1} with covariance ^T. 

T/ien for every ui,...,Ui G R the random processes {x n (uj,-),j = 
1, . . . , /} weakly converge in C([0; 1],R ) to i/ie l-point motion of the Harris 
flow with the local characteristic T. 

Proof. It follows from Lemma Q] and the invariance principle, that for ev- 
ery j = 1, . . . , / x n (uj,-) weakly converges in C([0; 1]) to the Brownian 
motion which starts from Wj. Then, it remains to prove that any limit 
point of {x n (uj, •), j = 1, . . . , /} coincides with the /-point motion of the 
Harris flow. Without loss of generality suppose that the whole sequence 
{x n (uj, •), j = 1, . . . , /} weakly converges. For a function / G C 3 (R) with 
bounded derivatives, consider the random processes 

Vn(t) = X n (Uj+i,t) ~ X n (Uj,t), 

z n {t) = f{y n {t)) - f(u j+ i - uj) - / (1 - T{y n (s)))f"{s)ds. 

Jo 

Following the known procedure (see for example [1]), it is easy to verify, 
that {z n ;n > 1} weakly converges to a certain martingale. Consequently 
the weak limit of y n satisfies the martingale problem for the operator 

Af(x) = (1 - T(x))^/(x). 

Since the martingale problem now has a unique solution [3] , then the weak 
limit of y n is the solution to the following Cauchy problem 

(dy(t) = ^2-2T(y(t))dw(t), 
\y(0) = u j+ i - uj. 
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The solution to this SDE has the strong Markov property. Consequently 
y is nonnegative for Uj + \ — Uj > 0. Hence, the weak limit of {x n (uj, •); j = 
1, . . . , 1} preserves the order. It remains to check the form of the joint cha- 
racteristic, which can be done in a standard way. The theorem is proved. 

□ 

The previous result is based on the uniqueness of a solution to SDE 
related to a stochastic flow. Now we consider the convergence of difference 
approximations to the n-point motions of the Arratia flow. Let us recall 
that Arratia's flow [3] is the Harris flow with the local characteristic T = 
U{o}. In this flow any two trajectories coalesce into a single one in finite 
time. 

Theorem 2. Suppose, that for every m > 1 x m is built upon a sequence 
n > 1} where independent identically distributed processes £™ have the 
covariance function T m which satisfies the Lipschitz condition. Define for 
m > 1 

2 - 2F m (x) 



C m = sup 



x 2 



1) lim^oc = 0, 

2) for every 5 > sup R \[_^ |r m (x)| ->• 0, m ->• oo, 

then the random processes {x m (ui, •),..., x m (ui, •); m > 1} weakly con- 
verge to the l-point motion of Arratia's flow starting from ui, . . . , u\. 

Proof. As in the proof of Theorem Q] we have the weak compactness of 
{(x m (ui, •),..., x m (ui, •); m > 1)} in C([0; 1], M 1 ) and the weak convergence 
of x m (ui, •) to a Wiener process. Consequently, for any limit point of 
{(x m (ui, •),... ,x TO (w;, •); m > 1)} it is enough to check the mutual char- 
acteristics and the order preserving property. For Ui < Ui + \ the difference 
process y m (t) = x m (ui + i,t) — x m (ui + i,t) are equidistributed with the dif- 
ference approximation v m to the solution of the SDE 

dy m (t) = \/2- 2T m (y m (t))dw(t), 
y m (0) = u i+1 - u { . 



It is known |T], that 



Esu P (v m (t)-y m (t)) 2 <C CmeCn 



[0;i] 



m 
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Note, that y m is nonnegative. Consequently, for every r > 
Pjinf y m < — rj = Pjinf v m < — r j — > 0, m — > oo. 

Vf 1 l [0;l] m J 

Hence the weak limit of any subsequence of {y m ;m > 1} is nonnegative. 
The completion of the proof can be done exactly as in the previous the- 
orem using martingale approximation and the fact that any nonnegative 
martingale remains at zero after hitting zero. The theorem is proved. □ 

3. Convergence of random maps 

In this section we will consider convergence of {x n ;n > 1} as random 
maps to corresponding maps from a stochastic flow. Let us begin with the 
case of smooth T. Define the sequence 

<V» = + -Un + l(<W), (5) 

v m 

where {^ n ; n > 1} is a sequence of independent stationary centered Gaus- 
sian processes with covariance function T satisfying the inequality 

VmGI: 1 - T(u) < Cu 2 

with some constant C. Define the Harris flow x corresponding to V. Note, 
that now x has a modification x(u,t),u G R, t G [0; 1] continuous with 
respect to both variables. Really, using the martingale inequality one can 
get that 

E sup (x(u, s) — x(v, s)) 2 < 

se[0;t] 

< 2(u - vf + 2E f (2 - 2T(x(u, s) - x{v, s)))ds < 

Jo 

< 2(u — v) 2 + 4C / sup (x(u, r) — x(v, r)) 2 ds. 

J0 re[0;s] 

Consequently, for some c 

E sup (x(u, t) — x(v, t)) 2 < c(u — v) 2 . 

te[0;l] 

This inequality together with the Kolmogorov condition gives us the desired 
property. 

The next statement asserts the convergence of our approximations to a 
stochastic flow in the case of smooth T. 
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Theorem 3. The random functions {x m = x™; m > 1} converge in distri- 
bution in the space C([a; b}) to the random function x for arbitrary interval 
[a;b]. 

Proof. The convergence of finite-dimensional distributions was proved in 
Theorem [H It remains to check the weak compactness of {x m ; m > 1}. For 
arbitrary u,ueRwe have 

E{x™ +l {u)-x™ +l {v)) 2 = E{x™{u)-x™{v)?+-E{2-2T{x™{u)-x™{v))) < 

m 

< E(x™(u) - x™{v)f + 2Cm- l E{x™{u) - x™{v)f. 
Consequently, 

orj 

E(x m (u) - x m {v)) 2 <(u- v)\l + —) m < e 2C (u - vf. 

m 

The obtained estimation gives the desired weak compactness. The theorem 
is proved. □ 

To obtain approximation of Arratia's flow we need some additional re- 
sults about the convergence of smooth stochastic flows to Arratia's flow. 
Let us consider the following SDE with the space-time white noise (Wiener 
sheet) W 

dz(u,t)= / tp(z(u,t)-p)W(dp,dt), 

Jr (6) 

z(u, 0) = u, ueK, 

where <p G Co°(R) and J R cp 2 (u)du = 1 (see P, [7] about equations of type 
(Ej)). All what we need here is a statement, that under our condition on (p 
the unique strong solution to (ED exists and is the Harris flow corresponding 
to the local characteristic 

T(u) = / (p(—p)(f(u — p)dp. 
Jr 

It was proved in [8], that the n-point motions of solutions z £ to (ED which 
corresponds to (f £ with the property supp ip £ C [— e; s] converge in distri- 
bution to the n-point motions of the Arratia flow when e — > 0. 
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Consider discrete approximations of z. For every n > 1 define 
zS(u) = u, 

4 + M = z n k (u) + / " / V(z n k (u) - p)W(dp, dt), (7) 

n 

k = 0, . . . , n — 1. 

It can be easily checked that every z% has a continuous modification. The 
next theorem gives a speed of convergence of 2™ to z(-,l) in the space 
C([0;1]). Define 

L 2 = / ^'{pfdp. 
Jr 

Theorem 4. There exist such positive constants C, C", C"\ that for every 
n>l 



E\\z n n - z(; 1)|| < ^= exp{(C"L 2 + C"'L A )e AL2 + L 2 }(L 2 + 1). (8) 



where \\ ■ \\ is the uniform norm in C([0; 1]). 
Proof. Consider for k = 1 , . . . , n 



E (zj(0) - z (o, ) 2 = B (^(O) - z (o, ^) ) 2 + 

+E ( n f (^(4-i(0) - p) - p(*(0, s) - p)) 2 ^s < 
J^- Jr 

n 

< E U_M - z U ^-T) V + I?E j^Jz'Um - 2(0, s)) 2 d S = 

= £ (^(o)-,(o,^)) 2 (i + ^) + 



n \ n 



2n 2 
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Consequently, 



£«(0)-z(0,l)) a < Je i2 



Note, that under our conditions on if, random functions {z^} and z have 
continuous derivatives with respect to the spatial variable. 

Let us denote by and y these derivatives. Then for k = 1, . . . , n 



Vliu) = yUiu) 1 + 




and 



Hence, 



dy(u,t)=y(u,t) / ^(u, *) - p)W(dp, dt). 



k 



k-: 



u)-y[u,-) = yt-iiu) -y[u,-) + 



n 



n 




+ / / -P)~ V( u i s)<p'{z(u, s) - p)]W{dp, dt) = 

fe-i 



= Vk-i(u) - y I u, 



k-l 



+ 



n 



+ 

+y [u 




yk-i( u ) -y\u. 



k-l 



n 



k-l 
n 

k-l\ 



i<p'(4-i(u)-p)-<p' (z (u, 



k-l 



n 



~P + 



n 



\z lu, 



k 'zl 

n 
k-l 



P) ~ <p'(z(%s) -p) I + 



+<p'(z(u, s) - p) (y (u, ~ y( u ' s ) 



W(dp, ds). 



Then 



£ (yk(u) -y ( u ^^j = E (vk-i(u) - y (u, 



k-l 



n 



+ 



+4£ 




yk-i( u ) -y( u 



k-i 



n 



ip\z n k _ l {u)-p) 2 dpds+ 



+4E [ n [ y ( u ,^—^\ y{zl_ l {u)-p)-ip' (z (u,^ ~ p) dpds+ 



k- lV ( ( k-1 



+4E / y [ u, \(p'(z[u, )—p )-(f'(z(u, s)—p)) dpds+ 

J^lJr \ n ) V V n J J 

+4E [ [ <p'(z(u,s) - p) 2 ( y ( u,- — - ) —y{u, s) J dpds < 
Jh=k ./to \ \ n 



k-l\\ 2 4 



+AEy L 2 J^(z ~ z ( w ' s )j ds + 

+4L 2 £ /" fj/ («, ^— ^ - s)^| ds. 



Jk-i \ \ n 
Note, that the processes z(u,t),t G [0; 1] and 

ri(t) = f f ip'(z(u, s) - p)W(dp, ds),te [0; 1] 
Jo Jr 

are continuous martingales with the characteristics 

(z(u,.))(t)=t, mt) = LH. 

Consequently, z(u, •) and 77 are Wiener processes. It follows from this that 

y(u,t) = exp{rj(t) - ^L 2 }. 



Hence 

Ey ( u, J / ( z ( u, ) — z(u, s) ) ds = 



2 -fc 



r, ( k~l\ 2 „ f'n ( ( k~l\ , \ 2 , 

= £/W U, £/ / Z M, — Z(U, S) ) ds < 

\ n J Jk_i \ \ n J J 

1 L 2 
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E lL{ y ( uk ~^)~ y{u ' s) ) ds= 



k 



2 



fc-i \ / fc- i 

n \ n 



= E I / r)(p'(z(u, r) — p)W(dp, dr) I <is = 




= L £ / / rydrds < 

Ik 
1 



J fc-1 

n n 



Furthermore 



< 'C 2 — < ' 



n 

In the last inequality the martingale property of x n and x was used. 
Finally one can get 

E Ll(u) - y L,^X\ < 



< E (yU(u) ~ V U "A V f 1 + ^ + + 5( £2 + 1 ' e 



Consequently, 

E(yl(u) - y(u, l)) 2 < 

< ^ exp{(c 5 L 2 + c 6 L 4 )e 4L2 + L 2 }(L 2 + 1). 
n 

To obtain an estimation for the uniform norm \\z™ — z(-, 1)|| we proceed as 
follows 

E\\z» - *(., 1)|| < E\zW) - *(0, 1)| + E f \y n n {u) - y(u, l)\du < 

Jo 

< ^= exp{(c 8 L 2 + c 9 L 4 )e 4L2 + L 2 }(L 2 + 1). 
'n 



The theorem is proved. □ 
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Obtained estimation can be used to get the convergence of the difference 
approximation to Arratia's flow. We will establish this convergence using 
the Levy-Prokhorov distance. Let us recall its definition. 

Definition 2. [9]. For two nondecreasing cadlag functions /, g on [0; 1] the 
Levy-Prokhorov distance is 

p(f,g)=M{E>0: V«G[0;1]: 

f(u — e) — e < g(u) < f(u + e) + e 
g(u — e) — e < f(u) < g(u + e) + e]. 

It is well-known p5] that the convergence in this distance is equivalent 
to the convergence at every point of continuity of the limit function. Also 
note that 

p(f,g) > d(f,g), 
where d(f,g) is the Skorokhod distance between / and g [2]. 
Take a function ijj G Cq° with supp^ C [—1; 1] such, that 

/ ^ 2 (u)du = 1. 
Jr 

For arbitrary e > define 

r e( w ) = - / lfje(p)lfje(u + P)dp. 
£ JR 

The parameter £ here is associated with the smoothness of T £ . In order 
to approximate the Arratia flow we have to take e — > 0. For independent 
Gaussian processes > 1} with the covariance {T £n } let us construct 

the sequences 

X n k+1 ( U ) = X H k (u) + —t n (x n k (u)). 

The next theorem shows that can be used to approximate the Arratia 
flow. 

Theorem 5. Suppose that £ n — > 0, n — > oo, 

— = o(llnn), n — > oo. 

£ 

Then the random functions converge weakly in D([0; 1]) to the value of 
the Arratia flow x(-, 1). 
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Proof. Consider the sequence of SDE 

dz En (u,t) = / tf) £ (z £n (u,t) - p)W(dp,dt). 
Jr 

As it was mentioned in the beginning of this section, for every ui, . . . , u m G 
[0; 1] {z En (ui, 1), . . . , z £n (u mj 1)) weakly converge to (x(u h 1), . . . , x(u m , 1)). 
Hence [10], z £n (-, 1) weakly converge to x(-, 1) in the Levy-Prokhorov dis- 
tance. For every n > 1 the sequence x", . . . , x™ is equidistributed with 
the discrete approximations to z £n from Theorem HI Consequently, x 1 ^ is 
equidistributed with x n such that 

E\\x n - z e ,X; 1)|| < ^=exp{(C"Ll + CTlSjJ*. + L E 2 J(4, + 1). 

where L\ n = j- J R i/j'(p) 2 dp. Hence 

E\\x n - z £n (-, 1)|| ->• 0, n ->■ oo. 
Since for continuous functions /, ^ the Skorokhod distance 

d{f,g) < \\f-g\\, 

then weakly converges to x(-, 1) in D([0; 1]). The theorem is proved. □ 
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